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There is experimental evidence that VLF signals propagating from west to east suffer 
less attenuation than from east to west. Earlier work treating the case of nonreciprocal 
propagation along the magnetic equator is extended in latitude. The nonreciprocity shown 
by \\R\\ for highly oblique propagation along the magnetic equator persists when the re- 
flection point moves towards a magnetic pole, but at the pole itself \\R\\ is reciprocal. To 
a first approximation j_R± is reciprocal at all magnetic latitudes. The conversion co- 
efficients \\R±, ±R\\, are greater for east-to-west propagation than for propagation in the 
opposite direction, except at a magnetic pole where they are equal. 

1. Introduction 

There is now a considerable amount of experimental evidence [Crombie, 1958; Taylor 
19G0J showing that VLF propagation is nonreciprocal, since VLF signals received from a west- 
erly quarter appear to suffer less attenuation than those received from the east. A qualitative 
explanation of such behavior in terms of interactions between the longitudinal movements of 
ionospheric electrons in the plane of incidence, and the transverse component of the earth's 
magnetic field was given by Crombie [1958]. 

Subsequently, Barber and Crombie [1959] calculated the reflection coefficient of a sharply 
bounded ionosphere and showed that this explanation was correct for the special case of waves 
polarized in the plane of incidence, propagating along the magnetic equator. This special 
case simplified the calculations since the magnetic field is horizontal and transverse to the 
plane of incidence. 

It has also been shown [Crombie, 1960] that the nonreciprocity which could be expected 
from the lack of reciprocity exhibited by the reflection coefficient does occur when mode theory 
is used. Wait [1960, 1961] and Wait and Spies [1960], using some coefficients derived from 
the numerical data of Johler [1961], have shown that mode propagation depends on the direc- 
tion of propagation. More recently, Dobrott and Ishimaru [1961] have also concluded that 
VLF propagation along the magnetic equator may be nonreciprocal. 

It is the purpose of this paper to derive formulas for the reflection coefficient of a sharply 
bounded ionosphere in the rather more general case where plane waves are incident from the 
west, or from the east at any magnetic latitude. 

2. Introductory Theory 

The starting point of this work is a paper b}^ Yabroff [1957a and b] in which, following 
an outline by Bremmer [1949], he derives perfectly general relationships between the fields in 
a sharply bounded ionosphere, and obtains the appropriate boundary equations. Yabroff 
considered the coordinate system shown in figure 1. The vertical direction is x and the earth's 
magnetic field H e is contained in the xz plane. The angle of incidence of the waves is &•, and 
the magnetic field of the earth is inclined to the vertical at an angle &. The incident, reflected 
and transmitted waves are described in the coordinate system x',y',z' in which for a wave 
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Figure 1. Coordinate system. 

The positive values of A and a» are measured clockwise 
from the positive x and z axes respectively. 



propagating in the a ?: ,ft direction, z r is in the a*,ft+ 90° direction; y f is directed in the a*+90°, 
90° direction and x' is in the direction at, Pi. The lower boundary of the supposedly uniform 
ionosphere is the xy plane. If z is magnetic north, then y is magnetic east. 
The incident plane wave F t can be written as 

J^=Re ~\ F oi exp j- (ct — x cos ft—;?/ sin ft sin a t —z sin ft cos a t ) > 

in which a>= wave angular frequency, c'is'the velocity of light and t is time. Since the boundary 
conditions must be satisfied it follows that the wave in the ionosphere F m , arising from this 
incident wave, must be of the form 

F m =T{e *{ F om exp j- (ct—Dx—y sin ft sin a t —z sin ft cos a ( ) V 

D may be regarded as the cosine of the (complex) angle of refraction in the ionosphere. 

Maxwell's equations and the equations of motion of an electron in a magnetic field can be 
written as 

dH "> 






m 



dt 



-eE-mvy-to(VXH e ) 



(1) 



in which, —e,m are the charge and mass of an electron, 
V is the electron velocity, 

v is the collision frequency, 

e , Mo are the permittivity and permeability of free space. 
From these three equations the three components of H and Fmay be eliminated giving three 

equations relating the components of E. In matrix form these equations are : 
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wh ore 

s=(l—jZ)/X h=Y/X, h L =h cos ft, h T =h sin ft 

X=Ne 2 /m€ o) 2 = (normalized plasma frequency) 2 . 
Y=n Q eH e lum= normalized gyro frequency. 
Z=v/o)= normalized collision frequency, 
a^^sin ft cos a i} a r =sin ft sin «*, a=sin ft. 

In order that i? be finite the determinant of (2) must be zero. Thus 

6J? 4 +63D 3 +6 2 Z? 2 +6 1 Z?+6 =0 (3) 

where 

& 4 =s(s 2 -A 2 )-s 2 +A£ 

b s =2a L h L h T 

b 2 =2s[(l-a 2 )h 2 -(s-l)(s-a 2 s-l)}-h 2 T (l-a 2 L )-(2-a 2 )h 2 L 

b x = — 2 (1 — a 2 )a L h L h T 

b =(s-l){[(l-a 2 )s-l] 2 -(l-a 2 ) 2 h 2 T cos 2 «,} 

-(l-a 2 )[(l-a 2 ),s— lj(^ sin 2 a«+Ai). (3a) 

The four complex values of Z> obtained from this quartic represent the two upgoing and two 
downgoing characteristic waves which can propagate independently in the medium. Values of 
D with positive real and negative imaginary parts represent upgoing waves. 

2.1. Boundary Equations 

Since the total tangential fields E y , E z , II y , and IL must be continuous across the boundary, 
the following boundary equations are obtained: 

From E y : 

Eyfi cos at+Ezfi cos ft sin oL t -\-E y , r cos ai—E z , T cos ft sin a i =P\E z i J rP 2 E z2 (4a) 

From E z : 

—E y fi sin ai-\-E z ,i cos ^ cosai—E y , r sin ai—E z , r cos ft cosa^Z^i+Z?^ (4b) 

From H y : 

E yfi cos ft sin a.i—E z ,i cos ai—E y , r cos ft sin a.i—E z , r cos a t =SiE z i-\-S 2 E z2 (4c) 

From i7 2 : 

JEy/j cos ft cos cKf+^Bg^ sin oii—E v t r cos ft cos c^+E^ sin ai=T\E z i-\-T 2 E z2 (4d) 

In these equations 

P^ / ^=M- 1 {-[l-a-^][a i a r -^] + [a^ + ^][a^ + ^]} 

(5) 

(6) 
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and 



where M is given by 

M -[ 1 -° , -«-^S)][ 1 -° i - c '-^]-[ o ' -^][ < " D+ i^] (7> 

R = Vo H x /E z =a T -a L P (8a) 

S=r, a H v IE z =a L Q-D (8b) 

T=r, H z /E 2 =DP-a T Q. (8c) 
In which 



'a— VW«o- 
2.2. Reflection Coefficient 2? 
This can be written in matrix form as 



R-- 



{&} < 9 > 



The first subscript denotes whether the electric field is parallel (||) or pendicular (_L) to the 
plane of incidence, and the second subscript refers in the same way to the electric field of the 
reflected wave. The components of R are defined as follows : 

\ft\\=E z , r \E zH (10a). Jt^E^/Eyn. (10c) 

$L L =E ylr \E zH (10b). L R ± =E yfr /E yn . (lOd) 

The components of R are obtained by solution of the boundary eqs (4) for the appropriate 
pairs of electric fields. 

3. Application to Propagation Along a Line of Magnetic Latitude 

The work of Yabroff given above is perfectly general and expressions for the reflection 
coefficient could be derived for any direction of propagation, and for any orientation of the mag- 
netic field. The labor involved would be excessive, however, and the resulting expressions 
cumbersome. The more practical approach to this general case is to use a computer for the 
solution of the boundary eqs (4). This has been done by Johler and Walters [1960] for propa- 
gation in the magnetic meridian. Johler [1961] has also calculated numerically, reflection coef- 
ficients for propagation at angles of 0°, 60°, 120°, 180°, 240°, 300° to the magnetic meridian, 
omitting the directions 90° and 270°, with w^hich this paper is concerned. 

However, in the case of propagation along a line of magnetic latitude the expressions 
obtained for the reflection coefficient are relatively compact as will be shown below, although even 
in this case the algebra is very tedious. 

For magnetic east-to-west or west-to-east propagation c^=90° or 270°. Thus, a L =0. 
Making this substitution in eqs (4) gives for the boundary equations: 

Ey: E/ i C-E z f r C=P 1 E zl +P 2 E z2 (11a) 

Hz: E/t+E/r^TiEa+TiE* (lib) 

Ez: -E v \-E y \=E zl +E z2 (lie) 

Hy: E/ i C-E/ r C=S 1 E zl +S 2 E z2 (lid) 

where C= cos p t . 

(8b) and (8c) reduce to T=BP-a T Q (12) 

S=-D (13) 
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while (5), (6), and (7) reduce to 

P=P'M- 1 =M-{i[l-a-|E^]+^(^+^)}A < 14) 

Q=Q'M-i=M->{ -j (a T D-J^y-f (i-Z>>-^)}^ (15) 

M=[l-a>-£^] [l-^- ii A_]-[a^+^] (16) 

The quantities M, P, Q, S, T, carry the same subscripts as the D which they contain. 

3.1. Solution of Boundary Equations 

The boundary eqs (11) can be solved in the following way for the components of R. \\R\\ 
is taken as an example. From (11a) and (lib) 

ii«ii-ji, ,/*,, t T l C+P l +(E a /E tl )(T i C+P 2 y {U) 

From (lie) and (lid), on putting E v ' i =0 il follows that 

Ea/E a = -(<7+A)/(<7+A). (18) 

Substitution of (18) in (17) gives 

li Pi=A- 1 f(7V-/ > l )(r+A)-(T 2 ^'-/ J 2)C'+A)l, (19) 

A= (T l C+P l )(C+ A) - (T 2 C+P 2 ) (C+DJ. (20) 

Similarly 

ll S x =A- I [2tf(A-A)L (2D 

x«ii=A- , [2<7(riP.-2y > ,)] > (22) 

± B ± =A-»[(Ti<7+Pi)(C r -A)-(OT 2 +P 2 )(<7-A)]. (23) 

Substitution for T u P 1} 7'.. and P 2 in these equations gives, after much work, the following 
relatively compact expressions for \\R\\ and xP-x 

,,/?„= A- [[{i[(l-J) <?AA--D8«7+A+A)]+^ [tf((7+A+A)-AA] }[^E^ 

.^^^l^^^^]! (24) 

x B x =Af'[^j[z?U(7-A-A)-AA(l-|)c]-^[AA-0((7+A+A)]}[ i ^ES) 

-l](A+iW4{i(2DJ+^+ 2 ^][(^-l)(«+IWW+^^^]|. (25) 

in which 

A '= A /[^2 (A-A)]=/i[(i-i) CDJX+Bl (C+A+A)] 

+ /^, [C(C+A+A)+AA] }[~E§)-i] (A+A) 
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and D\=\—a 2 — 1/.§, the value of D 2 when there is no magnetic field present. 
Equation (22) reduces to 

■*--*- , A[( i -0 i *+ i £] a 

Equation (21) reduces to 

ll-Rx=Af 1 ■2CM 1 M 2 (f^y 



(27) 



(28) 



When M X M 2 is expanded, it is found to contain terms involving D\-\-D\ and D\D\. These are 
respectively equal to — b 2 /b and bjb in eq (3). Substitution of these in (28) gives, after con- 
sid e rable reduction, 



\\Rj L =—A 1 



2h L 
s 2 -h/ 



[(i-i)« + ^] 



(29) 



which is identical to _J2||- 



4. Discussion of Reflection Coefficient Formulas 

These formulas have been checked under various limiting conditions to ensure that they 
conform with other published results. 



4.1. No Magnetic Field 
If the magnetic field is zero then \\R\\ and ±R±_ reduce to 

_(D 2 +a 2 )C-D Q 



\\R\\ 



±Ri 



\D 2 +a 2 )C+D 
C-D 



l± ~~O+D 
while ||B JL = _ L jB|| = 0. These limiting forms were given by Yabroff [1959b]. 



4.2. Horizontal Magnetic Field 
If the magnetic field is horizontal, h L =0, h T =h, and (19) reduces to 



iiBii= 



H>-9+^; 



[«(.-!> 



iIt^i 



7-Qvi- 



a 



fh T& tU'2 



C+(jDi+^f^ 



(30) 



From figure 1, h T is negative and for west-to-east propagation a T is also negative. In Barber 
and Crombie's paper [1959] their term equivalent to h T a T was taken as negative for west-to- 
east propagation, because of the different coordinate systems used by them. Bearing this 
in mind it is possible to show that (30) is identical with their earlier result. 
If in (25) h L =Q then (25) reduces again to 



±Rl~ 



C-D 1= C-Do 
'~C+D~C+D 



confirming the result of Barber and Crombie [1959]. It is easily seen from (27) and (29) that 
when h L =0 both j_JR|| and \\R± become zero. 
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4.3. Vertical Incidence, Vertical Magnetic Field 



Here h T , a T , a become zero and 0=1, 
Substitution in (24), (25), (27), (29), give 



DiA-1 n -n , A- A 



»Bii=j.«j.= (Z>i _ | . 1) {D2+1 y ±R\\=\\R±=j- (Di+lXD2+1) 
which are results previously obtained by Bremmer [1949], 

4.4. Reciprocity 

The expressions (24), (25), (27), and (29) all contain h T a T . This is positive for west-to- 
east propagation and negative for east-to-west propagation. Thus, provided that h T is not 
zero, all tin* reflection coefficient terms may exhibit nonreciprocity. 

5. Numerical Work 

In order to evaluate the expressions for the reflection coefficients it is necessary first of 
all to find D 1 and D 2 from (3). For propagation along a line of magnetic latitude a L =0 and 
(3) becomes a quadratic in D 2 . Thus 

bJD i +b 2 D 2 +b =0 i 

and 



^ 2 =[-&2±V^-46 6 4 ](26 4 )- 1 . 
Substituting a L =0 into eqs (3a) gives 

6 2 =2* [(l-a 2 )h 2 -«(s-l) Dl]-h 2 T -(2-a 2 )hl 
bl-4b b 4 =h 4 -2(l-anhVil+4s(s-l)Dlhl+(l-a 2 ) 2 hi 
h=s(« 2 -h 2 )-(s 2 -hi). 

The values of D u D 2 having positive real and negative imaginary parts are chosen since these 
represent the upgoing "ordinary" and "extraordinary" waves in the ionosphere. 
Some special cases are of interest in obtaining the results of section 4. 

For h=0; D 2 =D 2 =D 2 =l-a 2 -l/s. 

For h L =0; h T =h and D 2 =D 2 , D 2 =l-a 2 ' 



For h T =0; h L =h and a=Q, D 2 h2 =l- 



s 2 -h 2 -s 
1 



s±h 



Some approximate calculations have been made for daytime conditions to show the effect of 
the angle of dip of the earth's magnetic field on the four reflection coefficient terms using the 
formulas (24), (25), (27), and (29) developed above. The calculations are based on the fol- 
lowing parameters: 

o>o=10 6 corresponding to an electron density ~ 3 10/cc, 

v=10 7 /sec, co=10 5 /sec, a,«78°80', (7-0.2. 

The earth's magnetic field has been taken as 0.5 gauss and the inclination of the field to the 
vertical /3 e = 90°, 60°, 30°, and 0°. The results of these calculations are shown in table 1* 
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Table 1 





IP 


B|| 


± E 


JL 


\\ R ±- 


..1*11* 






West to 
east 


East to 
west 


West to 
east 


East to 
west 


West to 
east 


East to 
west 




90° 

60° 


0. 72/175° 
. 69/176° 
. 67/179° 
. 59/180° 


0. 33/175° 

. 38/176° 
. 42/183° 
.59/1; so 


0. 76/171° 
. 76/162° 
. 68/162° 
. 66/150° 


0. 76/171° 
. 76/162° 
. 68/160° 
. 66/150° 


0.0 

. 042/136° 
. 075/145° 
. 083/153° 


0.0 

. 081/192° 
. 122/177° 
. 083/153° 


Horizontal field. 


30°.... 




0° 


Vertical field. 



*For propagation in the Northern Hemisphere the phase angle of ii/2j_ and j_7?m should be reduced by 180°. 

It is clearly seen from this table that ||B||, \\Rj_ and j_R\\ all show nonreciprocal properties, 
whereas ±R± does not. In the case of \\R\\ the difference between east-to-west and west-to-east 
propagation is greatest at the magnetic equator where the magnetic field is horizontal, and 
gradually decreases as the reflection point moves nearer a magnetic pole. Nevertheless, the 
difference is very appreciable even when the angle of dip is 60° (£ c =30 o ). At the pole, where 
the magnetic field is vertical the two reflection coefficients are equal, since there is then no 
transverse field. 

At the magnetic equator \\Rj_ and ±R\\ are zero. When the magnetic field is inclined \\R± and 
±R\\ become finite, but again their values for east-to-west and west-to-east propagation differ. 
This difference again disappears at the magnetic pole. 

Despite the fact that eq (25) for ±R± contains a T h T , it seems from table 1 that, for the values 
of the parameters used, ±R± is essentially reciprocal. Thus in (25) the variations due to the 
sign of a T h T must, to a first approximation cancel. This presumably will not necessarily occur 
for other values of the parameters. 

The above numerical calculations have been made for highly oblique incidence, since long 
distance VLF propagation is of interest at present. However, the formulas are applicable at 
any angle of incidence. Thus the predictions of the formulas could be tested by making obser- 
vations of \\R\\ and/or ||JBj_ at points situated at equal distances to the west and east of a VLF 
transmitter. Alternatively, observation of the ground-interference pattern to the west and 
east of the transmitter could be made with airborne equipment. 

6. Conclusions 

This paper extends the earlier work of Barber and Crombie [1959] to the case of west-to- 
east and east-to-west propagation at any magnetic latitude. The nonreciprocity shown by 
l|i?n for highly oblique propagation along the magnetic equator persists when the reflection 
point moves towards a magnetic pole, but at the pole itself ||B|| is reciprocal. To a first approx- 
imation at least, ±R± is reciprocal at all magnetic latitudes. However, the conversion coef- 
ficients \\R±j ±R\\ y which are shown to be equal, are greater for east-to-west propagation than for 
propagation in the opposite direction, except at a magnetic pole where they are equal. 



The work described in this paper was done while the writer was on a short visit to the 
Boulder Laboratories of the National Bureau of Standards. The writer wishes to express his 
appreciation for this opportunity, particularly to Dr. J. R. Wait through whom it w T as arranged. 
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